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I present the status of the elements and parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
and summarise the related theoretical progress since Lepton-Photon 2003. One finds \V US \ = 0.2227 db 
0.0017 from K and r decays and \V c b\ = (41.6 db 0.5) • 10 -3 from inclusive semileptonic B decays. The 
unitarity triangle can now be determined from tree-level quantities alone and the result agrees well with 
the global fit including flavour-changing neutral current (FCNC) processes, which are sensitive to new 
physics. From the global fit one finds the three CKM angles 12 = 12.9° ±0.1°, <9 2 3 = 2.38° ±0.03° and 
013 = 0.223° ±0.007° in the standard PDG convention. The CP phase equals S 13 ~ 7 = (58.8 +|-|)° at 
la CL and 7 = (58.8 + is'^) a ^ 2cr CL. A major progress are first results from fully unquenched lattice 
QCD computations for the hadronic quantities entering the UT fit. I further present the calculation 
of three-loop QCD corrections to the charm contribution in — > tt+vV decays, which removes the 
last relevant theoretical uncertainty from the K — > irvV system. Finally I discuss mixing-induced CP 
asymmetries in b — > sqq penguin decays, whose naive average is below its Standard Model value by 
3a. 



1 Flavour in the Standard Model 

In the Standard Model transitions between 
quarks of different generations originate from 
the Yukawa couplings of the Higgs field to 
quarks. The non-zero vacuum expectation 
value v of the Higgs field leads to quark mass 
matrices M u and M d for the up-type and 
down- type quarks, respectively. The trans- 
formation to the physical mass eigenstate ba- 
sis, in which the mass matrices are diagonal, 
involves unitary rotations in flavour space. 
The rotation of the left-handed down-type 
quarks relative to the left-handed up-type 
quarks is the physical Cabibbo-Kobayashi- 
Maskawa (CKM) matrix V. It appears in the 
couplings of the W boson to quarks and is the 
only source of transitions between quarks of 
different generations. V contains one physi- 
cal complex phase, which is the only source of 
CP violation in flavour- changing transitions. 

Flavour physics first aims at the pre- 
cise determination of CKM elements and 
quark masses, which are fundamental pa- 
rameters of the Standard Model. The sec- 
ond target is the search for new physics, 



pursued by confronting high precision data 
with the predictions of the Standard Model 
and its extensions. To this end it is use- 
ful to distinguish between charged-current 
weak decays and flavour-changing neutral 
current (FCNC) processes. The determina- 
tion of CKM elements from the tree-level 
charged-current weak decays, discussed in 
Sect. 2, is practically unaffected by possible 
new physics. a By contrast, FCNC processes 
are very sensitive to virtual effects from new 
particles with masses at and above the elec- 
troweak scale, even beyond 100 TeV in cer- 
tain models of new physics. FCNC processes 
are discussed in Sect. 3. 

V can be parameterised in terms of three 
mixing angles #12, #23, #13 and one complex 
phase £13, which violates CP . Adopting the 
PDG convention 1 , in which V u d, V US: V c b and 
Vtb are real and positive, these parameters 



a Still new physics can be revealed if the 3x3 CKM 
matrix V is found to violate unitarity: One may then 
infer the existence of new (for example iso-vector) 
quarks which mix with the known six quarks. Further 
leptonic decays of charged mesons are tree-level, but 
sensitive to effects from charged Higgs bosons. 



C=(0,0) B=(1,0) 

Figure 1. Unitarity triangle (UT). 



Figure 2. \V u j\ and \V c j\, j = d,s,b, are determined 
from semileptonic decays. 



iS 13 



can be determined through 

= sin 0i2 cos 0i3, V ub = sin#i 3 e~ 

V c b = sin #23 COS #13. 

The Wolfenstein parameterisation 2 
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is an expansion of V in terms of A ~ 0.22 
to order A 3 . It shows both the hierarchy of 
the CKM elements and their correlations, like 
\V US \ — \V c d\ and \V cb \ ~ \V ts \. The apex of 
the standard unitarity triangle (UT), which 
is shown in Fig. 1 is defined by 3 



p + vq 



V ub V ud 



V ub V ud 



x , T „ T , (3) 

V ch Vcd 

(p, 77) coincide with (p, 77) up to corrections 
of order A 2 . With Eq. (3) and 

A = sin #12, AX 2 



sin 023, 



(4) 



the Wolfenstein parameterisation can be 
made exact 3 , that is V can be expressed in 
terms of (A, A,p, 77) to any desired order in 
A. In the following I always use the PDG 
phase convention and the exact definitions in 
Eqs. (3) and (4), with one exception: I ig- 
nore the small phase of — V cd (see Ref. 1 ), so 
that I can identify argF n * 6 = £13 with 7 and 
argV^ with the angle f3 of the unitarity tri- 
angle. This approximation is correct to 0.1%. 

The numerical results presented in the 
following have been prepared with the help of 



the Heavy Flavor Averaging Group (HFAG) 4 
and the CKMfitter 5 and UTFit 6 groups. 
CKMfitter uses a Frequentist treatment of 
theoretical uncertainties, while UTFit pur- 
sues a Bayesian approach, using flat prob- 
ability distribution functions for theoretical 
uncertainties. 



2 CKM elements from tree-level 
decays 

The standard way to determine the magni- 
tudes of the elements of the first two rows of 
V uses semileptonic hadron decays, depicted 
in Fig. 2. From Eq. (2) one realizes that an 
accurate determination of V us or V ud deter- 
mines V c d and V cs as well. Therefore mea- 
surements of semileptonic c — > d and c — > s 
decays are usually viewed as test of the com- 
putation of the hadronic form factors enter- 
ing the decay amplitudes. Charm decays are 
covered by Iain Stewart. 7 

2.1 V ud 

V u d can be determined from superallowed 
(0 + — > + ) nuclear f3 decay and from the (3 
decays n — > piV^j) and tt~ — > 7r° iV^j). 
Since no other decay channels are open, 
the semileptonic decay rate can be accessed 
through lifetime measurements. All three 
methods involve the hadronic form factor of 
the vector current: 

(f\uj^d\i), 
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where (ij) = (0+, 0+), (n,p) or (tt^tt ). 
The neutron (3 decay further involves the 
form factor of the axial vector current: 

( / l^7/i 75^1 i ) 

The form factors parameterise the long- 
distance QCD effects, which bind the quarks 
into hadrons. The normalisation of the vec- 
tor current is fixed at the kinematic point 
of zero momentum transfer Pi — Pf in the 
limit m u = rrid of exact isospin symmetry. 
The Ademollo-Gatto theorem 9 assures that 
corrections are of second order in the sym- 
metry breaking parameter (nid — m n )/A^ a ^, 
where A^ad is the relevant hadronic scale. 
No such theorem protects the axial form fac- 
tor (p 1^7^75^! n ), but the corresponding pa- 
rameter Ga can be extracted from asym- 
metries in the Dalitz plot. Experimentally 
the highest precision in the determination 
of V u d is achieved in the nuclear (3 decay, 
but n — > piVi(^) starts to become competi- 
tive. However, there is currently a disturbing 
discrepancy in the measurement of the neu- 
tron lifetime among different experiments. 8 
From a theoretical point of view progress in 
n — > p£Vi(j) and, ultimately, in the pristine 
7T~ — > 7T° iV^j) decay are highly desirable 
to avoid the nuclear effects of + — > + tran- 
sitions. On the theory side QED radiative 
corrections must be included to match the ex- 
perimental accuracy, recently even dominant 
two- loop corrections to n — > p£T>t(^) have 
been calculated. 10 

The world average for V u d reads 11 : 



0.9738 ± 0.0005 



(5) 



2.2 V UI 



V us can be determined from Kaon and r de- 
cays. The most established method uses the 
so-called K£3 decays K° ir-£ + v e , K° 

The decay rates schematically read 



and K + — > 7r°/i + ^. 
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Figure 3. y u jf° T (0). The horizontal band is 
the range quoted in Eq. (6). Courtesy of Vincenzo 
Cirigliano. 11 
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The hadronic physics is contained in 

(TT-{ P7r )\sru\K\ PK )) = 

a k + _ /r-°(o) 



X SU(2) 



1, 
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0. 



and QED corrections are contained in 
A^. The Ademollo-Gatto theorem 9 en- 
sures ff % "(0) = 1 + 0((m s - m d ) 2 /A 2 ad ). 
jK TV (^o) — 1 can be calculated with the 
help of Chiral Perturbation Theory (xPT) 12 , 
which exploits the fact that the pseudoscalar 
mesons are Goldstone bosons of a dynami- 
cally broken chiral symmetry of QCD. xPT 
amounts to a systematic expansion mp/ Ahad, 
^/Ahad, ^/Ahad an d the electroweak cou- 
pling e. Here p and M denote meson 
momenta and masses and mi is the lep- 
ton mass. There has been a substantial 
progress in the calculation of both A^ 13 and 
jK *- ( )i4 gince Lpm Significant effects of 
(D(e 2 p 2 ) QED corrections on differential dis- 
tributions were found; they must be included 
in Monte Carlo simulations. The value for 
V us f+ n (0) extracted from various experi- 
ments is shown in Fig. 3. The world average 
reads: 11 



fK"7T y 



0.2175 ±0.0008. 



(6) 
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Combining the results from x^T at order p 6 
and quenched lattice gauge theory (new) to 14 



n 

one arrives at 

V us 



K U 7T~ 



0.972 ±0.012 



0.2238 ±0.0029 



(7) 



from K£S. 

V us can also be determined from the Kfi2 
decay K + — > ^+^(7). 15 The hadronic quan- 
tity entering this decay is the Kaon decay 
constant Fk- Uncertainties can be better 
controlled in the ratio Fk / F n and one con- 
siders 

r(K+ - _ 



r(7^ 
V 2 

v us 

/•2 



M + ^(7)) 



F 2 K M K 



m, 



VL Fl Ml 



1 - - (C„ - C K ) 

7T 



with QED corrections C n — Ck — 3.0 ± 1.5. 
Using the result 16 F K / F^ = 1.210 ± 0.004 ± 
0.013 computed by the MILC collaboration 
with 2+1 dynamical quarks, one finds (with 
V ud from Eq. (5)): 



V us = 0.2223 ± 0.0026 



(8) 



from the Kji2 decay. This is astonishingly 
precise and K + — > /i + 1^(7) may constrain 
mass and couplings of a charged Higgs boson, 
which can mediate this decay as well. 15 

The third possibility to measure V us used 
hadronic r decays to the inclusive final state 
with strangeness |5| = 1. The experimental 
inputs are the ratios 

T AS=1 (r hadrons v T (j)) 



R TS — 



Rrd — 



r(r 
r AS=0 (r- 



ev e v T {i)) 
hadrons ^(7)) 



oc V 



oc V 



ud 



T(r -> ev e v T (^)) 

Here S is the strangeness. The optical the- 
orem allows to relate R TS ,d to the QCD 
current-current correlators II ^ d and 11^ d : 
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Figure 4. The optical theorem relates r AS_0 ' 1 (r — ► 
hadrons ^-(7)) to nj'j". The blob denotes the 

hadronic states contributing to nj'j". The leading 
term in the OPE is obtained by replacing the blob 
by a (u, d) or (u, s) quark loop and gluons to the de- 
sired order in a s . 



with z = s/M 2 = (p T -p VT ) 2 /M 2 . This re- 

T L 

lationship is depicted in Fig. 4. II s ' d can be 
computed through an operator product ex- 
pansion (OPE). The leading term is mass- 
less perturbative QCD, subleading operators 
entering IlJ' L are m 2 and m s {qq). The 
OPE amounts to an expansion in Aqcd/M t , 
m s /M T and a s (m T ). In the limit m s = of 
exact SU(3)f symmetry the ratio R TS /R T d 
would directly determine V 2 s /V 2 d . Hence it 
suffices to compute the (small) SU(3)f break- 
ing quantity 17 



5R r 



Rrd 

v u 2 d 



R TS 



With SRr = 0.218 ± 0.026 18 and experi- 
mental data from OPAL 19 one finds R T d = 
3.469 ± 0.014, R rs = 0.1694 ± 0.0049 and 
finally: 18 



V u . 



R T s 



V Rrd/\V ud \ 2 ~ SR r 

= 0.2219 ± 0.0033 exp ± 0.0009 th 
= 0.2219 ±0.0034. (9) 

The dominant source of uncertainty in 5R r , 
which enters Eq. (9) as a small correction, 
is from m s . In the near future it should be 
possible to improve on V us with data from 
BaBar and BELLE. 

In summary one finds an excellent consis- 
tency of the three numbers for V us from K£3, 
Kji2 and r decays. This is remarkable, since 
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= c X b 



Figure 5. OPE for B — ► X c ivg. The leading operator 
bb has dimension 3. 



the three methods use very different theoret- 
ical tools to address the strong interaction: 
Chiral perturbation theory, lattice gauge the- 
ory and the operator product expansion. The 
result nicely reflects the tremendous progress 
of our understanding of QCD at low energies. 
Averaging the results of Eqs. (7), (8) and (9) 
one finds: 



V us = 0.2227 ±0.0017 



(10) 



With V u d in Eq. (5) one can perform the first- 
row unitarity check 

vl ± v 2 d ± \v ub \ 2 - 1 ~ v 2 s ± v 2 d - 1 

= -0.0021 ±0.0012 



The Cabibbo matrix is unitary at the 1. So- 
level, just as at LP'03: 20 



V z +V 

v US 1 1 



ud 



-0.0031 ±0.0017 



2.3 V cb 



V c b can be determined from inclusive or ex- 
clusive b — > civ i decays. Exclusive decays are 
not discussed here. The analysis of the inclu- 
sive decay employs an OPE 21 , similarly to the 
determination of V us from r decay discussed 
in Sect. 2.2. The optical theorem relates the 
inclusive decay rate B — > X c £v£ to the imagi- 
nary part of the B meson self energy, depicted 
on the LHS of Fig. 5. The OPE matches 
the self energy diagram to matrix element 
of effective operators, whose coefficients con- 
tain the short-distance information associ- 
ated with the scale m& and can be calcu- 
lated perturbatively. Increasing dimensions 
of the operators on the RHS of Fig. 5 corre- 
spond to decreasing powers of m& in the co- 
efficient functions, so that the OPE amounts 
to a simultaneous expansion in Aqcd/^ and 



a s (m b ). Since (B\bb\B) = l+(D(A 2 QCD /m 2 b ) 
and there are no dimension-4 operators, non- 
perturbative parameters first occur at order 

a qgW™6- The y are 



a4 °c 



[B\bD 2 ± b\B] 



& oc (B\bia^G^b\B) 

[Iq, which parameterises the matrix element 
of the chromomagnetic operator, can be de- 
termined from spectroscopy. Hence to order 
A 2 Q CD /m 2 one only has to deal with the three 
quantities ra&, m c and /i 2 , which quantifies 
the Fermi motion of the b quark inside the B 
meson. 

The OPE can further be applied to cer- 
tain spectral moments of the B — > X£v£ de- 
cay, the distributions of the hadron invari- 
ant mass Mx and of the lepton energy. Fur- 
ther the same parameters govern different in- 
clusive decays, for instance also B — > X s j. 
Therefore there is a lot of redundancy in 
the determination of V c b, providing power- 
ful checks of the theoretical framework. The 
state of the art are fits to order Ag OT /m^, 
which involve 7 parameters. 22 The result of 
a global fit to hadron and lepton moments 
in B — > Xivi and photon energy moments 
in B X s <y from BaBar, BELLE, CDF, 
CLEO, DELPHI 23 can be seen in Fig. 6. It 
gives 

V c b = 41.6 ± 0.3 ex p ± 0.3oPE moments 

± 0.3oPEr s i 
= (41.6 ±0.5) • 10~ 3 (11) 

from inclusive B — > Xivg. 
24 \Vub\ 

I discuss the determination of \V u b\ from in- 
clusive B — > X u £v£ decays. Exclusive de- 
cays are discussed in. 7 In principle one could 
determine \V u b\ in the same way as V c b, if 
there were no background from B — > X c £vg 
decays. Its suppression forces us to impose 
cuts on the lepton energy Ei, the hadronic 
energy Ex, the hadron invariant mass Mx 
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Figure 6. Fit result for V c b vs. m^, which is defined 
in the kinetic scheme. 22 Fit and plot are courtesy of 
Oliver Buchmiiller and Henning Flacher. See also 23 . 



or a judiciously combination of them. Mx 
is too small for an OPE in the portion of 
phase space passing these cuts. Still some 
components of the hadron momentum Px are 
large. The description of inclusive B decays 
in this region involves the non-perturbative 
shape function 5, which is a parton distri- 
bution function of the B meson. At leading 
order in l/ra& the same S governs the pho- 
ton spectrum in B — > X s j and differential 
decay rates in B — > X u £p£. This allows us 
to extract S from B — > X s j for the use in 
B — > X u tp£. The goal to reduce the theoreti- 
cal uncertainty below 10% requires to under- 
stand corrections in both expansion parame- 
ters a s and AqcD/^b- E° r a correct treat- 
ment of radiative QCD corrections one must 
properly relate the differential decay rate dT 
to the shape function S. This is achieved 
by a factorisation formula, which has the 
schematic form: 24 

dT ocH I duo J (m b (P + - u)) S(u) 



Here H contains the hard QCD, associated 
with scales of order ra&. The jet function J 
and the shape function S contain the physics 



from scales of orders Mx ~ \/^b^QCD and 
A qcd, respectively. P + and P_ are defined 
as P± = E x T \Px V From A QCD < P+ ~ 
M x ~ ^rribAQCD < P- < rrib one re- 
alizes that one has to deal with a multi- 
scale problem, which is more complicated 
than B — > X c £v£. The second frontier of 
research in B — > X u tvg deals with sublead- 
ing shape functions Si, which occur at order 
l/rafo. They are different in B — > X u tvf> and 
P — > X s 7, but their moments can be related 
to OPE parameters like /i 2 , which gives some 
guidance to model these functions. 25 Mean- 
while an event generator for B — » X u tp£ de- 
cays is available, 26 with formulae which con- 
tain all available theoretical information and 
smoothly interpolate between the shape func- 
tion and OPE regions. It is pointed out that 
a cut on the variable P + , which is directly 
related to the photon energy in B — > X s 7, 
makes the most efficient use of the S(uj) ex- 
tracted from the radiative decay. 27 ' 26 Alter- 
natively one can eliminate S(uj) altogether by 
forming proper weighted ratios of the end- 
point photon and lepton spectra in B — > X s j 
and B — > X u tvg, respectively. 28 Using also 
the information from B — > XJLvg on and 
/i 2 the data from CLEO 29 , BELLE 30 and 
BaBar 31 combine to the world average 4 

V ub = (4.39 ± 0.20 exp ± 0.27 th , mb5/i 2 ) • 1(T 3 
= (4.39 ±0.34) • 10" 3 (12) 

from inclusive B — > X u tvg. Eq. (3) implies 
that |Kib/Vcb| defines a circle in the (p, 77) 
plane which is centered around (0,0). With 
Eqs. (11) and (12) its radius is constrained to 



0.45 ±0.04. (13) 



2.5 argKb 

7 = arg V* b can be determined from exclusive 
B — » { D q) X decays, where X denotes one or 
several charmless mesons. This method ex- 
ploits the interference of the tree-level b — » 
cuq and b — > ucq amplitudes, where q = d or 
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Figure 7. The Gronau-London-Wyler method com- 
( — ) 

bines the rates of B ± D° fi]^ for different 
final states fi. 



q = s. The prototype is the Gronau-London- 
Wyler (GLW) method 32 shown in Fig. 7. The 
decays B L>°X and B -> L>°X interfere, if 
both subsequent decays D° — > / and D° — > / 
are allowed. One needs four measurements 
to solve for the magnitudes of the b — > c 
and 6 — > u, their relative strong phase and 
their relative weak phase, which is the de- 
sired UT angle 7. For example one can com- 
bine the information of the branching frac- 
tions of B+ { D^\-^ K^^K* and B ± 
D° [— > 7r + 7r _ ]i ; C ± . This works with untagged 
non-flavour-specific decays as well: 33 E.g. the 
final state D° <j) does not reveal whether the 
decaying meson was a B s or B s . Still, when 
at least three pairs of B s — > D° [— > and 
£? s — > D° [— > branching fractions are 

measured, where f i = CPfi (and the f^s are 
not CP eigenstates), one has enough infor- 
mation to solve for 7. Since no flavour tag- 
ging is involved, the Tevatron experiments 
may contribute to these class of 7 determi- 
nations. The described determination of 7 
from tree-tree interference is modular, that 
is measurements in different decay modes can 
be combined, as they partly involve the same 
hadronic parameters. One should further 
first average the branching ratios from dif- 
ferent experiments and then determine 7 in- 
stead of averaging the inferred values of 7 ob- 
tained from different experiments. Combin- 
ing (almost) all B+ D°K+W data gives 



7 



(14) 



and the second solution 7 — 180° ~ —110°. 
This is 7 = arg V* b determined from the tree- 
level b — > ucs amplitude. 

Within the Standard Model b — > uud de- 
cays of tagged B° mesons are used to deter- 
mine the UT angle a. b — > uud decays involve 
both a tree and a penguin amplitude. The 
penguin component can be eliminated, if sev- 
eral decay modes related by isospin are com- 
bined, as in the Gronau-London method 34 
which uses B + — > 7r + 7r°, B° — > 7r + 7r~ and 
B° 7r°7r°. The B pn and 5 decay 
modes are better suited for the determina- 
tion of a, because the penguin amplitude is 
smaller. A combined analysis of the 7T7t, pir 
and pp systems gives 



a, 



exp 



(15) 



and the second solution a exp — 180° ~ —81°. 
The experimental result a exp could differ 
from the true a = Mg(-V* b Vtd/(V* b V ud )), if 
new physics alters the Bd — Bd mixing am- 
plitude. However, the influence from new 
physics is fully correlated in a exp and the 
CP asymmetry measured in b — > ccs decays. 
From the latter (see Eq. (21) below) we in- 
fer the Bd— Bd mixing phase 2/3 exp = (43.7 ± 
2.4)°. The Bd— Bd mixing phase cancels from 



the combination 27 
so that one obtains 



360° - 2a 



exp 



2A 



7 = (59ti 2 )° 



exp? 



(16) 



and the second solution 7 — 180° ~ —121°. 
Since the isospin analysis eliminates the pen- 
guin component, this is 7 = arg]/ n * 6 deter- 
mined from the tree- level b — > uud amplitude. 

The results in Eqs. (14) and (16) are in 
good agreement. Their naive average is 



7 



(63^)°. 



(17) 



The successful determination of a CP phase 
from a tree-level amplitude is a true novel re- 
sult compared to LP 1 OS. For the first time we 
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Figure 8. UT from tree quantities alone. The annulus 
is the constraint in Eq. (13) derived from \V u b\- The 
dark shadings correspond to 7 from Eqs. (14) and 
(16). Courtesy of Maurizio Pierini. 



can determine the UT from tree-level quanti- 
ties alone, the result is shown in Fig. 8. This 
is important, because the tree-level UT can 
only be mildly affected by new physics and 
therefore likely determines the true values of 
p and rj. 



q' u,c,t q 

Figure 9. Meson-antimeson mixing. (q,q f ) = (s,d), 
(6, d) and (6, s) for K— K mixing, B^— mixing and 
B s — B s mixing, respectively. 



are absent in generic extensions of the Stan- 
dard Model (like its supersymmetric gen- 
eralisations) making FCNC highly sensitive 
to new physics, probing scales in the range 
of 200 GeV to 100 TeV, depending on the 
model considered. This feature is a ma- 
jor motivation for the currently performed 
high-statistics experiments in flavour physics. 
Comparing different constraints on the UT 
from FCNCs processes and the tree-level con- 
straints discussed in Sect. 2 therefore pro- 
vides a very powerful test of the Standard 
Model. 



3 CKM elements from FCNC 
processes 

In the Standard Model FCNC processes are 
suppressed by several effects: First they only 
proceed through electroweak loops. Sec- 
ond they come with small CKM factors like 
\V ta \~ 0.04 and \V td \ ~ 0.01. Loops with an 
internal charm quark are further suppressed 
by a factor of m^/M^ from the GIM mecha- 
nism. Radiative and leptonic decays further 
suffer from an additional helicity suppres- 
sion, because only left-handed quarks cou- 
ple to W bosons and undergo FCNC tran- 
sitions. All these suppression mechanism are 
accidental, resulting from the particle content 
of the Standard Model and the unexplained 
smallness of most Yukawa couplings. They 



3.1 Meson-antimeson mixing 

K — K mixing, Bd — Bd mixing and B s — B s 
mixing are all induced by box diagrams, de- 
picted in Fig. 9. Each meson-antimeson sys- 
tem involves two mass eigenstates, their mass 
difference Am measures the magnitude of the 
box diagram and therefore constrains magni- 
tudes of CKM elements. The phase of box di- 
agram and thereby the phases of the CKM el- 
ements involved are constrained through CP- 
violating quantities. Tab. 1 shows the rela- 
tionship of the measurements to the CKM 
phenomenology. The quantities in the first 
two columns of Tab. 1 are well- measured and 
there is a lower bound on Am# s . 
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K — K mixing 


Bd — Bd mixing 


B s 


— B s mixing 


CP-conserving quantity: 


Attik 


Am Bd 




Ara# s 


CKM information: 


\v cs v cd \ 2 


\v tb v td \ 2 




\v tb v ts \ 2 


UT constraint: 


none 


R t = ^(l-p)*+rf* 




none 


CP-violating quantity: 


£k 


a™f(B d -> Jl^Ks) 


..mix 


(B s -> J/V>4>) 


CKM information: 


^{VtsV t * d ) 2 


sin(2/3) 




sin(2/3 s ) 


UT constraint: 


fj[(l — ~p) + const.] 


rj 
1-p 







Table 1. Relationship of meson-antimeson mixing to CKM and UT parameters. (3 = arg V^ d is one of the UT 
angles in Fig. 1 and (3 S = arg(— Vt s ) — \ 2 rj. 



3.2 e K 

While ex, which quantifies indirect CP vi- 
olation in K — > 7T7T decays, is measured at 
the percent level, its relationship to Im V t *j* oc 
r](l — p) is clouded by hadronic uncertainties 
in the matrix element 

( K° \ds v _ A ds v _ A \K° ) = % -f K M 2 K B K . 

This defines the hadronic parameter Bk-> 
which must be computed by non-perturbative 
methods like lattice QCD. Mk and fx are 
the well-known mass and decay constant of 
the Kaon. This field has experienced a major 
breakthrough since since meanwhile 

fully unquenched computations with 2+1 dy- 
namical staggered quarks are available. Us- 
ing MILC configurations the HPQCD collab- 
oration reports a new result 356 

B K (n = 2GeV) = 

0.618 ± 0.018 sta t ± 0.019 chira i 

extrapolation 
± 0.030discret. ± 0.130 p ert. matching 

= 0.618 ±0.136 (18) 

in the MS-NDR scheme. The conventionally 
used renormalisation scale and scheme inde- 
pendent parameter reads 

B K = 0.83 ±0.18 (19) 

The uncertainty from the perturbative 
lattice-continuum matching dominates over 

b In my talk I reported the preliminary value 
B K = 0.630 ± 0.018 s tat ±0.015 ch . extr . ± 0.030 disc . db 

0.130 p . match.- 



the statistical error and the errors from chiral 
extrapolation and discretisation in Eq. (18). 
This matching calculation was performed 
in 36 . The error in Eq. (18) is a conservative 
estimate of the unknown two-loop contribu- 
tions to this matching. If one instead takes 
twice the square of the one-loop result of 36 
as an estimate of the uncertainty, one finds 
0.036 instead of 0.130 in Eq. (18) and 

B K (fi = 2 GeV) = 0.618 ± 0.054 

B K = 0.83 ±0.07 (20) 

€k fixes 77(1—7?), so that it defines a hyperbola 
in the (p, 77) plane. 

3-3 V t d from Bd — Ba mixing 

The Bd — Bd mixing mixing amplitude in- 
volves the hadronic matrix element 

( B° \bdv-Abd V - A \B° ) = | M% d f Bd B Bd . 

Since the decay constant fs d is not mea- 
sured, the whole combination f% d BB d must 
be obtained from lattice QCD. The hadronic 
matrix element, however, cancels from the 
"gold-plated" mixing induced CP asymme- 
try aQp(Bd — > J/ipKs), which determines 
f3 = argV^ essentially without hadronic un- 
certainties. Combining all data from b — > ccs 
modes results in 37 ' 4 

sin(2/3) = 0.69 ± 0.03, cos(2/3) > 
=> wrg(±V t * d ) = /3 = (21.8 ±1.2)°. (21) 

The precisely measured Am# d = 0.509 ± 
0.004ps _1 is proportional to \V t d\ 2 f% d B Bd - 
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The HPQCD collaboration has computed 
f Bd = 216 ± 22MeV with 2+1 dynamical 


quantity 


central ± k_jL — lg 




staggered quarks. 38 This measurement is dis- 
cussed in detail in 7 . Combining this with B Bd 
from older quenched calculations results in 


P 
rj 


9 04 + - 035 

U.ZU^_ 

q oo 6 +0.021 
U.OOU _ .021 


+0.095 
-0.069 

+0.045 
-0.060 


fB d ^B Bd = (246 ±27)MeV, where B Bd = 


a (deg) 


98.4 If 1 , 


+ 16.8 
-11.8 


l.b2B Bd (n = m&) is the conventionally used 
scale and scheme independent variant of B Bd . 
Then from Am Bd alone we find 


P (deg) 
7 (deg) 


99 77 +0-87 
zz - ' ' -0.83 

58.8 ±|;1 


+ 1.92 
-2.04 

+ 11.2 
-15.4 


\v ub \ [io- 3 ] 


3.9o±8:g 


+0.29 
-0.24 


\V td \ = 0.0072 + 0.0008, 


\v td \ [io- 3 ] 


Q QQ +0.32 

°-^°-0.44 


+0.56 
-1.29 



where the error is reduced by a factor of 
2/3 compared to the old determination from 
quenched lattice QCD. 



3- 4 I Vtd \/\Vts\ from B — B mixing 

A measurement of the ratio ArriB d / Am Bs 
will determine |V^|/|V^| via 



Vtd 

Vts 



'Am Bd M Bs 



with 



the 



Am Bs V M B d 



hadronic 



quan- 



tity i = fB s yS B ^J(fB d yS Bd ) which equals 
£ = 1 in the limit of exact SU(3)f- A new un- 
quenched HPQCD result for f B Jf Bd 38 pre- 
sented in 7 can be used to refine the prediction 
for £. The lower bound Am Bs > 14.5 ps _1 
implies | Vt^/Vts | < 0.235 which constrains 
one side of the unitarity triangle: 



Rt = \/(l-p) 2 + rj 2 



V td 



VsX 



< 1.06 



3.5 Global fit to the unitarity triangle 

The result of a global fit of (p, rj) to state-of- 
the-art summer-2005 data is shown in Fig. 10. 
It uses B K = 0.85 + 0.02 + 0.07 where the first 
error is Gaussian and the second is scanned 
over according to the standard CKMfitter 
method 5 . For the remaining input see 5 . 
The fit output is summarised in this table: 



The output of the global fit agrees well 
with the pure tree-level determinations in 
Eqs. (12) and (17) and Fig. 8. 

We can use Eq. (1) to determine #13 = 
0.223° ± 0.007° from the fitted \V ub \ in the 
table. With Eq. (1) one finds 12 = 12.9° ± 
0.1° from Eq. (10) and <9 23 = 2.38° ± 0.03° 
from Eq. (11). Since 1 — cos #13 is negligibly 
small, the Wolfenstein parameters A and AX 2 
defined in Eq. (4) are simply given by V us in 
Eq. (10) and V c b in Eq. (11), respectively. 

3.6 K^itvV 

The rare decays K + — > tt+vV and Kl — > 
7r°uV provide an excellent opportunity to de- 
termine the unitarity triangle from s — > d 
transitions. With planned dedicated experi- 
ments (p, 77) can be determined with a sim- 
ilar precision as today from b — > d and 
b — > u transitions at the B factories. This 
is a unique and very powerful probe of the 
CKM picture of FCNCs. Br(K L tt°vV) 
is proportional to rj 2 and dominated by the 
top contribution. The theoretical uncer- 
tainty of the next-to-leading order (NLO) 
prediction 40 is below 2%. Br(K + — > n+vv) 
defines an ellipse in the (p, rj) plane and has a 
sizeable charm contribution, which inflicts a 
larger theoretical uncertainty on the next-to- 
leading order (NLO) prediction 41 , leading to 
(9(5 — 10%) uncertainties in extracted CKM 
parameters. Parametric uncertainties from 
V c b and m t largely drop out, if sin(2/3) is cal- 
culated from Br(KL — > 7r°uV) and Br(K + — > 
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Figure 10. LHS: The UT from a global fit to summer 2005 data. RHS: The first UT fit using theoretical 
expressions with NLO QCD corrections, performed in 1995. 39 At that time only \V u b\, ex and ArriB d could 
be used. Region la corresponds to a scan over la ranges of the input parameters. 



tt^uV). Therefore the comparison of sin(2/3) 
determined in Eq. (21) from the B system 
with sin (2/?) inferred from K — > i\vv consti- 
tutes a pristine test of the Standard Model. 42 
The impact of a future 10% measurements of 
these rare decay modes on the UT is shown 
in Fig. 11. 

The charm contribution is expanded in 
two parameters: m 2 K /m 2 c and a s (m c ). The 
calculations of 0(m 2 K /m 2 c ) corrections was 
recently completed, finding a 7% increase 
of Br(K + — > tt+vV) with a small residual 
uncertainty. 43 A new result are the next-to- 
next-to-leading order (NNLO) QCD correc- 
tions to the charm contribution. 44 This three- 
loop calculation reduces the theoretical error 
from unknown higher-order terms well below 
the parametric uncertainty from m c . The 
branching ratio is now predicted as 

Br(K + n + Vv) = (8.0 ± 1.1) • 10" 11 . 

At NNLO one finds the following reduced 
theoretical uncertainties for parameters ex- 
tracted from Bv(Kl — > ifivv) and Br(K + — > 



^-^ = 0.010, 5sm(2(3) = 0.006, 57 = 1.2° 

\ V td\ 

4 CP violation in b — » s penguin 
decays 

Within the Standard Model the mixing- 
induced CP asymmetries in b — > sqq pen- 
guin amplitudes are proportional to sin(2/3) eff 
which equals sin(2/3) in Eq. (21) up to small 
corrections from a penguin loop with an up 
quark. In b — > suu decays there is also a 
color-suppressed tree amplitude. In any case 
the corrections are parametrically suppressed 
by \V ub V us /(V cb V cs )\ ~ 0.025. The experi- 
mental situation is shown in Fig. 12. A naive 
average of the measurements of Fig. 12 gives 

sin(2/?) eff = 0.51 ±0.06, 

which is below the value of sin(2/3) from tree- 
level b — > ccs decays in Eq. (21) by 3<r. More- 
over QCD factorisation finds a small and pos- 
itive correction to sin(2/3 e ff )— sin(2/3) from up- 




-1 -0.5 0.5 1 1.5 2 
P 



Figure 11. Unitarity triangle from a future measure- 
ment of Bt(Kl — > 7r°i/77) and Br(K + — > ir+vV) with 
10% accuracy. The current UT from the 2005 global 
fit 5 is overlaid. The comparison of this UT con- 
structed from s d FCNCs with the UT found from 
b decays probes the CKM origin of FCNCs precisely 
and in a unique way. 



sin(2p e )/sin(2(^i ) F™P 

PRELIMINARY 



Figure 12. sin(2/3) eff from various penguin decays. 
The small vertical yellow band is sin(2/3) from 
Eq. (21). 4 



quark effects. 45 While the significance of the 
deviation has decreased since the winter 2005 
conferences, the mixing-induced CP asymme- 
tries in b — > sqq decays stay interesting as 
they permit large effects from new physics. 
While in Bd decays the needed interference 
of a Bd and Bd decay to the same final state 
requires a neutral K meson in the final state, 
b — > sqq decays of B s mesons go to a flavour- 
less ssqq state, so that the desired CP effects 
can be studied in any final state. Hence B s 
physics has the potential to become the "El 
Dorado" of b — > sqq penguin physics. 
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DISCUSSION 

Luca Silvestrini (Rome and Munich): 

Maybe rather than saying that QCD 
cannot explain sin2/3 eff in b — > s pen- 
guins, one should say that a particu- 
lar model of power supressed corrections 
due to Beneke and Co. cannot do it, but 
this is not a model-independent state- 
ment. If you just want to use data, and 
you say you do not know anything about 
power corrections, I do not think that 
you can infere anything from that plot. 

Ulrich Nierste: The parametric suppres- 
sion of the up-quark pollution by 

\V ub V us /(V cb V cs )\ ~ 0.025 is undisputed. 
Further the leading term in the \jm\y ex- 
pansion of sin 2/? eff — sin(2/3) can be reli- 
ably computed and results in the finding 
of Ref. 45 that sin2/3 eff - sin(2/3) is small 
and positive for the measured modes. It 
is true that the size of the modeled power 
corrections is currently widely debated. 
Yet I am not aware of any possible dy- 
namical QCD effect in two-body B de- 
cays which is formally 0(l/mb), large in 
magnitude and further comes with the 
large strong phase needed to flip the sign 
of sin2/f ff -sin(2/3). 



